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II. A Memoir on the Theory of Matrices. By Arthur Catlby, Esq., F.B.S. 



Beceiyed December 10, 1857, — Eead January 14, 1858. 



The term matrix might be used in a more general sense, but in the present memoir I 
consider only square and rectangular matrices, and the term matrix used without quali- 
fication is to be understood as meaning a square matrix; in this restricted sense, a set 
of quantities arranged in the form of a square, e. g. 

( a , b , c ) 
a', b', d 
a", b", c" 

is said to be a matrix. The notion of such a matrix arises naturally from an abbreviated 
notation for a set of linear equations, viz. the equations 

~K=ax -\-by +cz, 

Y =a'x -\-b'y +dz, 

Z =a"x+b"y+c"z, 

may be more simply represented by 

(X, Y, Z)=(.a, b , c Jx,y,z), 
d, b', d 
a", b", c" 

and the consideration of such a system of equations leads to most of the fundamental 
notions in the theory of matrices. It will be. seen that matrices (attending only to those 
of the same order) comport themselves as single quantities ; they may be added, 
multiplied or compounded together, &c. : the law of the addition of matrices is pre- 
cisely similar to that for the addition of ordinary algebraical quantities; as regards their 
multiplication (or composition), there is the peculiarity that matrices are not in general 
convertible; it is nevertheless possible to form the powers (positive or negative, 
integral or fractional) of a matrix, and thence to arrive at the notion of a rational and 
integral function, or generally of any algebraical function, of a matrix. I obtain the 
remarkable theorem that any matrix whatever satisfies an algebraical equation of its 
own order, the coefficient of the highest power being unity, and those of the other 
powers functions of the terms of the matrix, the last coefficient being in fact the deter- 
minant ; the rule for the formation of this equation may be stated in the following con- 
densed form, which will be intelligible after a perusal of the memoir, viz. the determi- 
mdccclviii. D 



18 



ME. A. CATLEY ON THE THEOKY OF MATEICES. 



nant, formed out of the matrix diminished by the matrix considered as a single quantity 
involving the matrix unity, will be equal to zero. The theorem shows that every 
rational and integral function (or indeed every rational function) of a matrix may be 
considered as a rational and integral function, the degree of which is at most equal to 
that of the matrix, less unity ; it even shows that in a sense, the same is true with respect 
to any algebraical function whatever of a matrix. One of the applications of the 
theorem is the finding of the general expression of the matrices which are convertible 
with a given matrix. The theory of rectangular matrices appears much less important 
than that of square matrices, and I have not entered into it further than by showing 
how some of the notions applicable to these may be extended to rectangular matrices. 

1. For conciseness, the matrices written down at full length will in general be of the 
order 3, but it is to be understood that the definitions, reasonings, and conclusions 
apply to matrices of any degree whatever. And when two or more matrices are spoken 
of in connexion with each other, it is always implied (unless the contrary is expressed) 
that the matrices are of the same order. 

2. The notation 

( a , I , c Jx, y, z) 

a!, b', c' 

a", b", c" 

represents the set of linear functions 

((a, b, c\x, y, z), (a', b', c'\x, y, z), (a", b", c"Jx, y, z)), 

so that calling these (X, Y, Z), we have 

(X, Y, Z)=( a , b ,c Jx, y, z) 

«', b',c' 

a", b", c" 

and, as remarked above, this formula leads to most of the fundamental notions in the 
theory. 

3. The quantities (X, Y, Z) will be identically zero, if all the terms of the matrix 
are zero, and we may say that 

( 0, 0, ) 

0, 0, I 

0, 0, o I 

is the matrix zero. 

* 
Again, (X, Y, Z) will be identically equal to (x, y, z), if the matrix is 

( 1, 0, ) 
0, 1, 
0, 0, 1 
and this is said to be the matrix unity. We may of course, when for distinctness it is' 
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required, say, the matrix zero, or (as the case may be) the matrix unity of such an order^ 
The matrix zero may for the most part be represented simply by 0, and the matrix 
unity by 1. 

4. The equations 

(X, Y, Z)=( a,b,c X*> Vi *)> (*'. Y', Z')=(-« , d , y Xx,y,z) 



d, b', d 
a", W, c" 



«', /3', y' 
«", j3", 7" 



(X+X', Y+Y, Z+Z')=( a +« , J +0 , c +7 », y, «) 

a'+u', b'+fi', d -+-y' 

y+«", #"+/3", c "+y" 

and this leads to 

( a +a, 5 +/3,,C +y )=(a,i,c ) + ( a , /3 , y ) 



«'+«', J'+/3', c'+y' 
a" +a ", #' + |3" ? c" +y " 



a', J', c' 



«\ /3', y' 
«", /3", y" 



as the rule for the addition of matrices ; that for their subtraction is of course similar 
to it. 

5. A matrix is not altered by the addition or subtraction of the matrix zero, that is, 
we have M+0=M. 

The equation L=M, which expresses that the matrices L, M are equal, may also be 
written in the form L — M=0, i. e. the difference of two equal matrices is the matrix 
zero. 

6. The equation L=— M, written in the form L+M=0, expresses that the sum of 
the matrices L, M is equal to the matrix zero, the matrices so related are said to be 
opposite to each other; in other words, a matrix the terms of which are equal but 
opposite in sign to the terms of a given matrix, is said to be opposite to the given 
matrix. 

7. It is clear that we have L+M=M+L, that is, the operation of addition is com- 
mutative, and moreover that (L+M)+N=L+(M+N)==L+M+N, that is, the opera- 
tion of addition is also associative. 

8. The equation 
(X, Y, Z)=( a , b , c ~$mx, my, mz) 

a!, b', d J 

a", b", d' 
written under the forms 

(X, Y, Z)=m( a , b , c \w, y, z)=( ma , mb , mc \x, y, z) 



a', b', d 

a", b", d' 



ma', mb' , md 
ma", mb", md' 



d2 
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gives 



w( a , b , c )=( ma , m# , mc ) 



a', 5', d 
a", 5", c" 



ma', mb' , md 
ma", mb", md' 



as the rule for the multiplication of a matrix by a single quantity. The multiplier m 
may be written either before or after the matrix, and the operation is therefore com- 
mutative. We have it is clear m(L+M)=mL+wM, or the operation is distributive. 

9. The matrices L and mL may be said to be similar to each other ; in particular, if 
m=l, they are equal, and if m= — 1, they are opposite. 

10. We have, in particular, 

w( 1, 0, )=f m, 0, ), 
0, 1, 0, m, 

0, 0, 1 0, 0, m 

or replacing the matrix on the left-hand side by unity, we may write 

m=( m, 0, ), 
0, m, 
0, 0, m 

The matrix on the right-hand side is said to be the single quantity m considered as 
involving the matrix imity. 

11. The equations 

(X, Y, Z)=( a, b , c Jw,y, z), (oc, y, z)=( a. , /3 , y T§, n, Q, 



give 



a', b', d 
a", b", d' 

(X,Y,Z)=( A, B, C #, *,£)=( «> *> G 

a', b', d 
a", b", d' 



1 



*', P, y 7 
*", P\ y" 



«', /3', y' 
«", i3'', / 



A', B', C 
A", B", C" 
and thence, substituting for the matrix 

( A, B, C ) 

A', B', O 

A", B", C" 
its value, we obtain 

((*, J, <?X«, «',«"), («, J,o Xl3,/y,j3"), («,*,cXy,y',/))=( a »*' c I*>^»y ) 



(«", 5", c"X«, «', «"), («", «", C'XA P, P'), («", *", ^"Xy, •, y") 



a',b',d 
a", b", d' 



*',P,y< 
a", /3\ / 
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as the rule for the multiplication or composition of two matrices. It is to be 
observed, that the operation is not a commutative one; the component matrices may 
be distinguished as the first or further component matrix, and the second or nearer 
component matrix, and the rule of composition is as follows, viz. any line of the com- 
pound matrix is obtained by combining the corresponding line of the first or further 
component matrix successively with the several columns of the second or nearer com- 
pound matrix. 

12. A matrix compounded, either as first or second component matrix, with the matrix 
zero, gives the matrix zero. The case where any of the terms of the given matrix are 
infinite is of course excluded. 

13. A matrix is not altered by its composition, either as first or second component 
matrix, with the matrix unity. It is compounded either as first or second component 
matrix, with the single quantity m considered as involving the matrix unity, by multi- 
plication of all its terms by the quantity m : this is in fact the before-mentioned rule 
for the multiplication of a matrix by a single quantity, which rule is thus seen to be a 
particular case of that for the multiplication of two matrices. 

14. We may in like manner multiply or compound together three or more matrices : 
the order of arrangement of the factors is of course material, and we may distinguish 
them as the first or furthest, second, third, &c, and last or nearest component matrices, 
but any two consecutive factors may be compounded together and replaced by a single 
matrix, and so on until all the matrices are compounded together, the result being inde- 
pendent of the particular mode in which the composition is effected ; that is, we have 
L.MN=LM.N=LMN, LM.NP=L.MN.P, &c, or the operation of multiplication, 
although, as already remarked, not commutative, is associative. 

15. We thus arrive at the notion of a positive and integer power 1/ of a matrix L, 
and it is to be observed that the different powers of the same matrix are convertible. 
It is clear also thdXp and q being positive integers, we have I/.I/=I/ + ", which is the 
theorem of indices for positive integer powers of a matrix. 

16. The last-mentioned equation, I/.I/=I/ +? , assumed to be true for all values what- 
ever of the indices p and q, leads to the notion of the powers of a matrix for any form 
whatever of the index. In particular, I/.L°=I/ or L°=l, that is, the Oth power of a 
matrix is the matrix unity. And then putting jb=1, q= — 1, or p= — 1, q=l, we have 
L.L -1 =L~ l .L=l ; that is, IT 1 , or as it may be termed the inverse or reciprocal matrix, 
is a matrix which, compounded either as first or second component matrix with the 
original matrix, gives the matrix unity. 

17. We may arrive at the notion of the inverse or reciprocal matrix, directly from the 
equation 

(X, Y, Z)=( a , b , c Jx, y, z), 

a 1 , V, d 
a", b", d 1 
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in fact this equation gives 

(x, y, z)=( A, A', A" XX, Y, Z)=(( a , b , c )-*IX, Y, Z), 



B, B', B" 

C, C, C" 



«', 5', c' 
a!', b",d' 



and we have, for the determination of the coefficients of the inverse or reciprocal matrix, 
the equations 

( A, A', A" Y a , b , e )=( 1, 0, ), 



B, B', B" 

C, C, C" 



a', V, d 
a", b", d' 



0, 1, 
0, 0, 1 



( a , b , g 



X 



a', b', d 
a", b", c" 



A, A, A" )=( 1, 0, ), 



B, B', B" 0, 1, 

C, C, C" 0, 0, 1 

which are equivalent to each other, and either of them is by itself sufficient for the com- 
plete determination of the inverse or reciprocal matrix. It is well known that if V 
denote the determinant, that is, if 

a , b , c ; 
a!, b 1 , d 
a!', b", d' 
then the terms of the inverse or reciprocal matrix are given by the equations 

, &c. 



A=i 
V 


1, 0, 

0, b', d 


V 


0, 1, o 
a', 0, d 




0, b", d' 




a", 0, d' 



or what is the same thing, the inverse or reciprocal matrix is given by the equation 

( a, b, c )-> _I( B.V, B„,V, B a „V ) 



a', b', d 
a", b\ d' 



d A V, B 4 ,V, 3 4 „V 
Z C V, B C ,V, o,,V 



where of course the differentiations must in every case be performed as if the terms a, b, 
&c. were all of them independent arbitrary quantities. 

18. The formula shows, what is indeed clear a priori, that the notion of the inverse 
or reciprocal matrix fails altogether when the determinant vanishes : the matrix is in 
this case said to be indeterminate, and it must be understood that in the absence of 
express mention, the particular case in question is frequently excluded from considera- 
tion. It may be added that the matrix zero is indeterminate ; and that the product of 
two matrices may be zero, without either of the factors being zero, if only the matrices 
are one or both of them indeterminate. 
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19. The notion of the inverse or reciprocal matrix once established, the other nega- 
tive integer powers of the original matrix are positive integer powers of, the inverse or 
reciprocal matrix, and the theory of such negative integer powers may be taken to be 
known. The further discussion of the fractional powers of a matrix will be resumed in 
the sequel. 

20. The positive integer power D" of the matrix L may of course be multiplied by 
any matrix of the same degree, such multiplier, however, is not in general convertible 
with L ; and to preserve as far as possible the analogy with ordinary algebraical func- 
tions, we may restrict the attention to the case where the multiplier is a single quan- 
tity, and such convertibility consequently exists. We have in this manner a matrix cL m , 
and by the addition of any number of such terms we obtain a rational and integral func- 
tion of the matrix L. 

21. The general theorem before referred to will be best understood by a complete 
development of a particular case. Imagine a matrix 

M=( a, b ), 

I c, d 
and form the determinant 

a— M, b 

c, d—M 

the developed expression of this determinant is 

W-(a+d)W + (ad-bc)W; 
the values of M 2 , M 1 , M° are 

( a*+bc , b(a+d) ), ( a, b ), ( 1, ), 



c(a-\-d), d 2 -\-bc 



c, d 



0, 1 



and substituting these values the determinant becomes equal to the matrix zero, viz. we 
have 



a— M, b 
c , d 



•M 



= ( a'+bc , b(a+d) )-(a+d) ( «, b )+( a d-bc) ( 1, ) 



c(a-\-d), d?+bc 



c, d 



that is, 



I 0, 1 I 

=( (cf+bc)— (a+d)a+{ad— be), b(a+d)—(a+d)b ) = ( 0, ), 

I c(a-\-d)—(a-\-d)c , d*+bc—(a-\-d)d-\-ad—bc \ | 0, j 

:0 



a— M, b 
c , d— M 
where the matrix of the determinant is 

( a, b )-M( 1, ), 

I £? d I I 0, 1 J 

that is, it is the original matrix, diminished by the same matrix considered as a single 
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quantity involving the matrix unity. And this is the general theorem, viz. the deter- 
minant, having for its matrix a given matrix less the same matrix considered as a single 
quantity involving the matrix unity, is equal to zero. 

22. The following symbolical representation of the theorem is, I think, worth no- 
ticing: let the matrix M, considered as a single quantity, be represented by M, then 
writing 1 to denote the matrix unity, M.l will represent the matrix M, considered as a 
single quantity involving the matrix unity. Upon the like principles of notation, 1 . M 
will represent, or may be considered as representing, simply the matrix M, and the 
theorem is 

Det. (T.M-M.1)=0. 

2 S. I have verified the theorem, in the next simplest case, of a matrix of the order 3, 
viz. if M be such a matrix, suppose 

M=( «, b, c ), 

d, e,f 
g, h, i 
then the derived determinant vanishes, or we have 

a— M, b , c = 0, 

d , e-U, f 

9 , h , «— M 
or expanding, 

M 3 — (a+e-\-i)W-\-(ei+ia+ae— fh— eg— bd)M— (aei+bfg+cdh— afh— bdi—ceg)=0 ; 
but I have not thought it necessary to undertake the labour of a formal proof of the 
theorem in the general case of a matrix of any degree. 

24. If we attend only to the general form of the result, we see that any matrix what- 
ever satisfies an algebraical equation of its own order, which is in many cases the mate- 
rial part of the theorem. 

25. It follows at once that every rational and integral function, or indeed every 
rational function of a matrix, can be expressed as a rational and integral function of an 
order at most equal to that of the matrix, less unity. But it is important to consider 
how far or in what sense the like theorem is true with respect to irrational functions of 
a matrix. If we had only the equation satisfied by the matrix itself, such extension 
could not be made; but we have besides the equation of the same order satisfied by the 
irrational function of the matrix, and by means of these two equations, and the equa- 
tion by which the irrational function of the matrix is determined, we may express the 
irrational function as a rational and integral function of the matrix, of an order equal at 
most to that of the matrix, less unity; such expression will however involve the coeffi- 
cients of the equation satisfied by the irrational function which are functions (in number 
equal to the order of the matrix) of the coefficients assumed unknown, of the irrational 
function itself. The transformation is nevertheless an important one, as reducing the 
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number of unknown quantities from n* (if n be the order of the matrix) down to n. To 
complete the solution, it is necessary to compare the value obtained as above, with the 
assumed value of the irrational function, which will lead to equations for the determina- 
tion of the n unknown quantities. 

26. As an illustration, consider the given matrix 

M=( a, b ) 

c, d I 

and let it be required to find the matrix L=\/M. In this case M satisfies the equation 

M 2 — (a+d)M+ad—bc=0 ; 
and in like manner if 

L=( «, (3 ) 

I 7-> & I 
then L satisfies the equation 

L 2 -(a+5)L+a&-£y=0; 

and from these two equations, and the rationalized equation L 2 =M, it should be possible 
to express L in the form of a linear function of M : in fact, putting in the last equation 
for L 2 its value (=M), we find at once 

L =^[M+(«$-/3y)], 

which is the required expression, involving as it should do the coefficients oc +5, aft— fly 
of the equation in L. There is no difficulty in completing the solution ; write for short- 
ness a+&=X, «(5— j3y=Y, then we have 

L=( «, (3 ) = ( a+Y 4 ), 
' XX 



5 



c 



d+Y 
XX 



and consequently forming the values of a+l and a&— /3y, 

a+d+2Y 



X= 



Y= 



X ' 

{a + Y){d+Y)-bc 



X 2 

and putting also a+d=P, ad—bc=Q, we find without difficulty 

X= x /p+2 N /Q, 
Y =x /Q, 

and the values of a, |3, y, & are consequently known. The sign of \/Q is the same in 
both formulae, and there are consequently in all four solutions, that is, the radical -y/M 
has four values. 

MDCCCLVIII. e 
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27. To illustrate this further, suppose that instead of M we have the matrix 
M 2 =( a, b )•=( <f+bc , b(a+d) ), 



c, d 



c[a-\-d), d 2 -\-bc 



so that L 2 =M 2 , we find 

P=(a+dy-2(ad-bc), 

Q=(ad—bcf, 
and thence \/Q= +(ad— be). Taking the positive sign, we have 

Y=ad— be, 

X=±(a+d), 
and these values give simply 

L= + ( a, b )=±M, 

\ c, d \ 
But taking the negative sign, 

Y = — ad-\-bc, 

X=±*/(a-d) 2 +4bc, 

and retaining X to denote this radical, we find 

L==( a s -ad+2bc b(a + d) ), 
X ' X 



e(a+d) d^—ad+2bc 
X ' X 



which may also be written 



L _«+<* (a, b ) 2(ad-bc) ( 1, ), 



c, d 



X 



0, 1 



or, what is the same thing, 



T « + J 2(ad— be) . 

L -~x~ M x , 



and it is easy to verify a posteriori that this value in fact gives L S =M 2 . It may be 
remarked that if 

M 2 =( 1, o y=i, 

0, 1 | 

the last-mentioned formula fails, for we have X=0; it will be seen presently that the 
equation L 2 =l admits of other solutions besides L=+L The example shows how the 
values of the fractional powers of a matrix are to be investigated. 

28. There is an apparent difficulty connected with the equation satisfied by a matrix, 
which it is proper to explain. Suppose, as before, 

M=( c, b ), 

I c, d \ 
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so that M satisfies the equation 



a— M, b 

c , d— M 



= 0, 



or 

M 2 -(a+d)M+ad-fo=0, 

and let X,, X ;/ be the single quantities, roots of the equation 



•X, b 

, d-X. 



=0 



or 

X 2 -(«+^)X+«^-^=0. 

The equation satisfied by the matrix may be written 

(M-X ; )(M-X„)=0, 

in which X ; , X„ are to be considered as respectively involving the matrix unity, and it 
would at first sight seem that we ought to have orie of the simple factors equal to zero, 
which is obviously not the case, for such equation would signify that the perfectly inde- 
terminate matrix M was equal to a single quantity, considered as involving the matrix 
unity. The explanation is that each of the simple factors is an indeterminate matrix, 
in fact M— X, stands for the matrix 

( a-X„ 6 ), 

I c , d—X t I 

and the determinant of this matrix is equal to zero. The product of the two factors is 
thus equal to zero without either of the factors being equal to zero. 

29. A matrix satisfies, we have seen, an equation of its own order, involving the 
coefficients of the matrix ; assume that the matrix is to be determined to satisfy some 
other equation, the coefficients of which are given single quantities. It would at first 
sight appear that we, might eliminate the matrix between the two equations, and thus 
obtain an equation which would be the only condition to be satisfied by the terms 
of the matrix ; this is obviously wrong, for more conditions must be requisite, and we 
see that if we were then to proceed to complete the solution by finding the value of the 
matrix common to the two equations, we should find the matrix equal in every case 
to a single quantity considered as involving the matrix unity, which it is clear ought 
not to be the case. The explanation is similar to that of the difficulty before adverted 
to, the equations may contain one, and only one, common factor, and may be both of 
them satisfied, and yet the common factor may not vanish. The necessary condition 
seems to be, that the one equation should be a factor of the other ; in the case where 
the assumed equation is of an order equal or superior to the matrix, then if this equation 
contain as a factor the equation which is always satisfied by the matrix, the assumed 
equation will be satisfied identically, and the condition is sufficient as well as necessary : 

e2 
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in the other ease, where the assumed equation is of an order inferior to that of the 
matrix, the condition is necessary, but it is not sufficient. 

30. The equation satisfied by the matrix may be of the form M M =1 ; the matrix is 
in this case said to be periodic of the nth. order. The preceding considerations apply to 
the theory of periodic matrices ; thus, for instance, suppose it is required to find a matrix 
of the order 2, which is periodic of the second order. Writing 

M=( a, b ), 

c, d I 
we have 

W-(a+d)M+ad-bc=0, 

and the assumed equation is 

M 2 -1 = 0. 

These equations will be identical if 

a-\-d=0, ad—bc=—l, 

that is, these conditions being satisfied, the equation M 2 — 1 = required to be satisfied, 
will be identical with the equation which is always satisfied, and will therefore itself 
be satisfied. And in like manner the matrix M of the order 2 will satisfy the condition 
M 3 — 1=0, or will be periodic of the third order, if only M 3 — 1 contains as a factor 

W-(a+d)M+ad-bc, 
and so on. 

31. But suppose it is required to find a matrix of the order 3, 

M=( a, b, c ) 
d, e,f 

g, h, i 
which shall be periodic of the second order. Writing for shortness 

«-M, b , c = -(M 3 -AM 2 +BM-C), 

d , e-M, f 

g , h i—U 

the matrix here satisfies 

M 3 -AM 2 +BM-C=0, 

and, as before, the assumed equation is M 2 — 1=0. Here, if we have 1+B=0, A-fC=0, 
the left-hand side will contain the factor (M 2 — 1), and the equation will take the form 
(M 2 — 1)(M+C)=0, and we should have then M 2 — 1=0, provided M+C were not an 
indeterminate matrix. But M+C denotes the matrix 

( a+C, b , c ) 

d , e+C, f 

g , h , «*+c 

the determinant of which is C 3 +AC 2 +BC+C, which is equal to zero in virtue of 
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the equations l-f-B=0, A-fC=0, and we cannot, therefore, from the equation 
(M 2 — 1)(M+C)=0, deduce the equation M 2 — 1=0. This is as it should be, for the 
two conditions are not sufficient, in fact the equation 

M 3 =( a'+bd+cg, ab+be+ch, ac+bf+ci )=1 

da-{-ed-\-fg, db-\-e 2 +/A, dc-\-ef-\-fi 

ga-\-hd-\-ig, gb-\-he-{-ih, gc-\-hf-\-i 2 

gives nine equations, which are however satisfied by the following values, involving in 
reality four arbitrary coefficients ; viz. the value of the matrix is 

* =( _^__ . Zl + 1% " (/3+7) ; 



) 



a + /3 + y 

'— (y+«)f t, ~ 



P 



•(*+.)£ 



(« + /3)j*v- 



a 7 

'■ + P+y a+#+y ' « + /3 + y 

so that there are in all four relations (and not only two) between the coefficients of the 
matrix. 

32. Instead of the equation M M — 1=0, which belongs to a periodic matrix, it is in 
many cases more convenient, and it is much the same thing to consider an equation 
~M. n —k=0, where k is a single quantity. The matrix may in this case be said to be 
periodic to a factor pres. 

33. Two matrices L, M are convertible when LM=ML. If the matrix M is given, 
this equality affords a set of linear equations between the coefficients of L equal in 
number to these coefficients, but these equations cannot be all independent, for it is 
clear that if L be any rational and integral function of M (the coefficients being single 
quantities), then L will be convertible with M ; or what is apparently (but only appa- 
rently) more general, if L be any algebraical function whatever of M (the coefficients 
being always single quantities), then L will be convertible with M. But whatever the 
form of the function is, it may be reduced to a rational and integral function of an order 
equal to that of M, less unity, and we have thus the general expression for the matrices 
convertible with a given matrix, viz. any such matrix is a rational and integral function 
(the coefficients being single quantities) of the given matrix, the order being that of the 
given matrix, less unity. In particular, the general form of the matrix L convertible 
with a given matrix M of the order 2, is L=«M+0, or what is the same thing, the 
matrices 

(M), ( < V ) 



will be convertible if a'—d': b' 



| c, d 
d=a — d 



d! 



b:c. 
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34. Two matrices L, M are skew convertible when LM=— ML; this is a relation 
much less important than ordinary convertibility, for it is to be noticed that we cannot 
in general find a matrix L skew convertible with a given matrix M. In fact, con- 
sidering M as given, the equality affords a set of linear equations between the coeffi- 
cients of L equal in number to these coefficients ; and in this case the equations are 
independent, and we may ehminate all the coefficients of L, and we thus arrive at a 
relation which must be satisfied by the coefficients of the given matrix M. Thus, sup- 
pose the matrices 

(■«,*)>( < V ) 



c, d 



d, d' 



are skew convertible, we have 

( a, b ¥ a', if )=( aa'+bd, ab'+bd' ), 



c, d d, d' 



ca'+dd, cV +dd! 



( a', V )( a, b )=( aa'+b'c, a'b+b'd ), 



d, d' c, d 



da+d'c, c'b-\-d'd 



and the conditions of skew convertibility are 

2aa'+bd+b'c =0 
b'(a+d)+b(a'+d')=0 
d(a+d)+c(a'+d')—0 
2dd'+bd+b'c =0 

Eliminating a', b', d, d', the relation between a, b, c, d is 



2a, 


c , 


h , 




b, 


a-\-d, 


• 


b 


c, 


• . 


a-\-d, 


c 




c , 


b , 


2d 



=0. 



which is 

(a+df(ad—bc)=0. 

Excluding from consideration the case ad— bc=0, which would imply that the matrix 
was indeterminate, we have a-\-d=0. The resulting system of conditions then is 

a+d=0, a'+d'=0, aa'+bd -\-b'c+dd'=0, 

the first two of which imply that the matrices are respectively periodic of the second 
order to a factor pres. 

35. It may be noticed that if the compound matrices LM and ML are similar, they 
are either equal or else opposite ; that is, the matrices L, M are either convertible or 
skew convertible. 
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( a, b ), ( a, c ) 



c, d 



b, d 



are said to be formed one from the other by transposition, and this may be denoted by 
the symbol tr. ; thus we may write 

( a, c )=tr.( a, b ), 



b, d 



c, d 



The effect of two successive transpositions is of course to reproduce the original matrix. 

37. It is easy to see that if M be any matrix, then 

(tr.M) p =tr.(M*), 
and in particular, 

(tr.M)~ 1 =tr.(M- 1 ). 

38. If L, M be any two matrices, 

tr.(LM)=tr.M. tr.L, 
and similarly for three or more matrices, L, M, N, &c, 

tr. (LMN)= tr. N. tr. M. tr. L, &c. 

40. A matrix such as 

( a, h, g ) 

h, b,f 

*/> c 
which is not altered by transposition, is said to be symmetrical. 

41. A matrix such as 

( 0, „, -it, ) 

-f, 0, X 

p, —x, o 

which by transposition is changed into its opposite, is said to be skew symmetrical. 

42. It is easy to see that any matrix whatever may be expressed as the sum of a sym- 
metrical matrix, and a skew symmetrical matrix; thus the form 

( « , A+», g— V* ) 
h-v, b ,/+X 

9+t*> /- x ' € 

which may obviously represent any matrix whatever of the order 3, is the sum of the 
two matrices last before mentioned. 

43. The following formulae, although little more than examples of the composition of 
transposed matrices, may be noticed, viz. 

( a, I J (8, <! )=( # 2 +# 2 ? ac+bd ) 



c, d d, b 



ac-\-bd, c a +^ 2 
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which shows that a matrix compounded with the transposed matrix gives rise to a sym- 
metrical matrix. It does not however follow, nor is it the fact, that the matrix and 
transposed matrix are convertible. And also 

( a, c T a, b Y a, c )=( a 3 +bcd+a(b a +c i ), (f+aM+^+d 8 ) ) 



b, d 



c, d 



b, d 



b 3 +acd+b(a?+d*), d 3 +al>c+d(b i +<?) 



which is a remarkably symmetrical form. It is needless to proceed further, since it is 
clear that 

( a, c Y «, b Y a, c X a, b )=(( a, c ^ a, b )) 2 



b, d <?, d b, d c, d 



b, d c, d 



44. In all that precedes, the matrix of the order 2 has frequently been considered, 
but chiefly by way of illustration of the general theory ; but it is worth while to develope 
more particularly the theory of such matrix. I call to mind the fundamental properties 
which have been obtained, viz. it was shown that the matrix 

M=( a, b ), 

c, d 



W~(a+d)M+ad-bc=0, 
(a, b ), ( a', b< ), 



c, d 



V, d' 



satisfies the equation 
and that the two matrices 

will be convertible if 

a'—d ! ;b':c=a—d:b:c, 

and that they will be skew convertible if 

a+d=0, a'+d'=Q, aa'+bd+b'c+dd'^O, 

the first two of these equations being the conditions in order that the two matrices may 
be respectively periodic of the second order to a factor pres. 

45. It may be noticed in passing, that if L, M are skew convertible matrices of the 
order 2, and if these matrices are also such that L 2 =— 1, M 2 =— 1, then putting 

N=LM=— ML, we obtain 

L 2 =-l, M 2 =-l, N»=-l, 

L=MN=-NM, M=NL=-NL, N=LM=-ML, 

which is a system of relations precisely similar to that in the theory of quaternions. 

46. The integer powers of the matrix 

M=( a, b ), 

c, d J 

are obtained with great facility from the quadratic equation; thus we have, attending 
first to the positive powers, 
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M*=z(a+d)M-(ad— be), 

W=[(a+dy-(ad-bcJ]M-(a+d)(ad-bc), 

&c, 

whence also the conditions in order that the matrix may be to a factor pres periodic of 

the orders 2, 3, &c. are 

a+d = 0, 

(a+d) 2 -(ad-bc)=0, 

Sec; 

and for the negative powers we have 

(ad—bc)M.~ l =—M.+(a+d), 

which is equivalent to the ordinary form 

(ad—bc)M- l =( d, -b); 

I — c, ® I 

and the other negative powers of M can then be obtained by successive multiplications 
with M". 

47. The expression for the nth. power is however most readily obtained by means of 
a particular algorithm for matrices of the order 2. 

Let h, b, c, J, q be any quantities, and write for shortness R= —b?—ibc ; suppose also 
that h', b', d, J', q 1 are any other quantities, such nevertheless that h' :b' : d=h:b: c, and 

write in like manner R'=— h' 2 — ib'd. Then observing that — t^> — t=» —7= are respect- 

h' V c' 
ively equal to -j=j t -^ -^==, the matrix 

( / h \ 2bJ 

75 ' K cotq+ ~h) 

contains only the quantities J, q, which are not the same in both systems ; and we may 
therefore represent this matrix by (J, q), and the corresponding matrix with h', b', d, J', q' 
by (J', q 1 ). The two matrices are at once seen to be convertible (the assumed relations 
h! : b' : d= h:b:c correspond in fact to the conditions, a'—d' : b' : d= a—d : b : c, of con- 
vertibility for the ordinary form), and the compound matrix is found to be 



ysin q sin q ' ■* ' * / 



And in like manner the several convertible matrices (J, q), (J', q 1 ), (J", q") &c. give 
the compound matrix 



/ 8in(g+g |+ / ) jjfJ „ \ 

\smqsmq' sin 9".. ' 1 ■ a ■ 1 j 



MDCCCLVIII. 
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48. The convertible matrices may be given in the first instance in the ordinary form, 
or we may take these matrices to be 

( a, b ), ( «', V ), ( «", b" ) &c. 



c, d I I </, dl 



c", d" 



where of course d—a : b : c=d'—a! : b' : d=d"—a" : b" : c"=&c. Here writing h=d—a, 

d+a 

~7r: 



and consequently R= — (d— of— ibc, and assuming also 1=^/11 and cot q=— 7==, and 



in like manner for the accented letters, the several matrices are respectively 
and the compound matrix is 

49. When the several matrices are each of them equal to 

( «, b )' 
I c, d I 
we have of course q=q l =:q".., R=R'=R".., and we find 



I ^5 ^ 



=(SHhmr.«): 



or substituting for the right-hand side, the matrix represented by this notation, and 
putting for greater simplicity 



S?f(i^)-=(K/S)L, „xL=S» (is / E) .-, 
we find 

( a, b ) n =(±l(s/R cot nq— (d— a)), U 



) 



c, d 



Lc ,%L(*S~R cot nq-\-(d— afj 



where it will be remembered that 

R=— (d— of— ibc and cotg , = —^■> 

the last of which equations may be replaced by 

d+a + */ — R 



cos 



q+*/—lsinq=- 



2 */ad—bc 



The formula in fact extends to negative or fractional values of the index n, and when n 
is a fraction, we must, as usual, in order to exhibit the formula in its proper generality, 
write q+Zrmr instead of q. In the particular case w=-J, it would be easy to show the 
identity of the value of the square root of the matrix with that before obtained by a 
different process. 
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50. The matrix will be to a factor pres, periodic of the nth order if only sin nq=0, 
that is, if q= — (to, must be prime to n, for if it were not, the o\ der of periodicity would be 

not n itself, but a submultiple of n); but cosgs= Q ■ , -^ , and the condition is therefore 

(d+af-4(ad-bc) cos 2 ^=0, 
or as this may also be written, 

(iP+ffl 2 — lad cos V-Abc cos 2 — =0, 

a result which agrees with those before obtained for the particular values 2 and 3 of 
the index of periodicity. 

51. I may remark that the last preceding investigations are intimately connected with 

dX ~f" b 

the investigations of Babbagb and others in relation to the function tpx— ^zj 
I conclude with some remarks upon rectangular matrices. 

52. A matrix such as 

( a, !>, c ) 
a', V, d | 
where the number of columns exceeds the number of lines, is said to be a broad matrix ; 
a matrix such as 

( a, b ) 

' a', V 
a", b" 
where the number of lines exceeds the number of columns, is said to be a deep matrix. 

53. The matrix zero subsists in the present theory, but not the matrix unity. 
Matrices may be added or subtracted when the number of the lines and the number of 
the columns of the one matrix are respectively equal to the number of the lines and the 
number of the columns of the other matrix, and under the like condition any number 
of matrices may be added together. Two matrices may be equal or opposite the one to 
the other. A matrix may be multiplied by a single quantity, giving rise to a matrix of 
the same form ; two matrices so related are similar to each other. 

54. The notion of composition applies to rectangular matrices, but it is necessary that 
the number of lines in the second or nearer component matrix should be equal to the 
number of columns in the first or further component matrix; the compound matrix will 
then have as many lines as the first or further component matrix, and as many columns 
as the second or nearer component matrix. 

55. As examples of the composition of rectangular matrices, we have 

( a, b, c % a 1 , b',d, a )=((«, b, cja', d, i'), (a, b, cjb',f,j'){a, b, cjd, g', K), (a, b, cjd', h\ I 1 ) ), 
d,e, f | | d , /', </, h' \(d,e, fXa\ d, i>), (d, e, fjb\ f\ J'M * fX<*> <?'> K % (*» *> fX*> h '> l ') \ 

f2 
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( a d X < V, o', «?))=( («, W. A (a, *X*V/)> («, W A («, W *) )• 



b e 
of 



d,f,g',h< 



(b, eXct, d), (b, ejb',f), (b, eX<?, 9% (*, «X*> h') 

56. In the particular case where the lines and columns of the one component matrix 
are respectively equal in number to the columns and lines of the other component 
matrix, the compound matrix is square, thus we have 

( a, b, c X a ', <$' ) = ( («» *» c X a '> h '> <0» («. 5 > C X^'> e '>/0 ) 



d, «, / 



and 






( «', eP X «. h c )=( 



5', e' 
o',f 



d, e, f 



(d, e, fXa', V, d), (d, e, /£<?, d, /') 



(a', ^'X«, d), (a\ d'Xb, e), (a 1 , d'Xc, f) ) 



(*', dXa, d), (V, dXh e), (b>, dXc, f) 
(d, fXa, d), (d, fXb, e), (d, fXc f) 

The two matrices in the case last considered admit of composition in the two different 
orders of arrangement, but as the resulting square matrices are not of the same order, 
the notion of the convertibility of two matrices does not apply even to the case in 
question. 

57. Since a rectangular matrix cannot be compounded with itself, the notions of 
the inverse or reciprocal matrix and of the powers of the matrix and the whole resulting 
theory of the functions of a matrix, do not apply to rectangular matrices. 

58. The notion of transposition and the symbol tr. apply to rectangular matrices, the 
effect of a transposition being to convert a broad matrix into a deep one and reciprocally. 
It may be noticed that the symbol tr. may be used for the purpose of expressing the 
law of composition of square or rectangular matrices. Thus treating (a, b, c) as a 
rectangular matrix, or representing it by (a, b, c), we have 






tr.( a', V, d )=( 



o, 



and thence 



( a, b, c )tr.( a', b', d )=(a, b, cX a' )=(«> b, cXd, V, d), 
, , , , ^ , v 

d 
so that the symbol 

(a, b, cX«', *\ d) 

would upon principle be replaced by 

( a, b, c ) tr. ( a', b', d ): 

ill I 
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it is however more convenient to retain the symbol 

(a, b, c£a', b\ d). 

Hence introducing the symbol tr. only on the left-hand sides, we have 

( a, b, c )tr.( a', if, d )=((«, b, cja!, b', d), (a, b, cjd!, d, /') ), 
\d,e,f\ | d>, d, f I I (d, e, fja', b\ d), (d, e, f\d>, d, f) \ 

or to take an example involving square matrices, 

( a, b )tr.r a', V )=( (a, b\a!, V), (a, b\d', d) ), 

I d, e I I d', d I | (d, dja\ b'), (d, e\d\ d) \ 

so that in the composition of matrices (square or rectangular), when the second or 
nearer component matrix is expressed as a matrix preceded by the symbol tr., any line 
of the compound matrix is obtained by compounding the corresponding line of the first 
or further component matrix successively with the several lines of the matrix which 
preceded by tr. gives the second or nearer component matrix. It is clear that the terms 
' symmetrical' and ' skew symmetrical' do not apply to rectangular matrices. 



